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Fine details of the excitation spectrum of the two-dimensional spin-gap magnet PHCC are revealed 
by electron spin resonance investigations. The values of anisotropy parameters and the orientations 
^ ' of the anisotropy axes are determined by accurate measurements of the angular, frequency-field and 

temperature dependences of the resonance absorption. The properties of a spin-gap magnet in the 
vicinity of critical field are discussed in terms of sublevel splittings and p-factor anisotropy. 
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I. INTRODUCTION 

Piperazinium hexachlorodicuprate (C4Hi2N2)(Cu2Cl6), abbreviated as PHCC, was studied recently as a good 
test example of a two-dimensional Heisenberg antiferromagnet.— ^— Inelastic neutron scattering experiments revealed 
that each magnetic ion interacts with at least six in-plane neighbours. The strongest exchange integral value is 
approximately 1.3 meV. The complex geometry of the exchange bonds can be envisioned as a set of strongly coupled 
antifcrromagnctic ladders running along the a axis (from here on we use the axes notation of Refi) of the triclinic 
crystal while some of the weaker exchange bonds are frustratedi^ The magnetic ground state of PHCC is a non- 
magnetic singlet separated from the excited triplet states by an energy gap of A = 1.02 meVji 

Due to the presence of the gap PHCC remains in the disordered spin-liquid state down to temperatures well below 
the exchange integral scale. At low temperatures the magnetic properties of this system can be considered as those 
of a gas of weakly-interacting triplets. The application of an external field lowers the energy of one of the triplet 
substates and, at a certain critical field, a quantum phase transition to a field-induced antifcrromagnctic state occurs. 
The closing of the gap and the formation of magnetic order above the critical field of ~ 8 T was confirmed by magnetic 
measurements as well as by calorimetric and neutron scattering techniques^. Additionally, this magnet was used as a 
test probe for the quasiparticle breakdown in a quantum spin liquid^. 

Magnetization measurements revealed an anisotropy of the critical field^. This anisotropy points to deviations from 
the Heisenberg model. Such deviations are important for the physics of a spin-gap magnet, especially around the 
quantum phase transition. Anisotropic interactions are known to cause a gap reopening in fields above the phase 
transition and they could also affect the quasiparticle decay by allowing otherwise forbidden processes. 

In this paper we present results of electron spin resonance (ESR) investigations of PHCC. Our study reveals the 
presence of anisotropic interactions leading to a g-factor anisotropy and a splitting in the triplet sublevels. We have 
determined magnitudes of these anisotropic effects and the orientations of the corresponding axes. Additionally, we 
address the question of the correct implementation of a macroscopic theorjs^i^ to the case of a spin-gap magnet with 
a relevant ^-factor anisotropy. 



II. EXPERIMENTAL DETAILS AND SAMPLES 

The samples were grown from saturated solutions of PHCC. The saturated solution was prepared by adding of 
piperazine C4H10N2 (from Sigma Aldrich), dissolved in minimal amount of concentrated HCl to CuCl2 •2H2O (99.99% 
purity from Sigma Aldrich), dissolved in minimal amount of concentrated HCl, at a 4:1 molar ratio of CuCl2 to 
C4H10N2. 

The main set of samples for the ESR measurements was grown by cooling the filtered saturated solution in the 
refrigerator. After several days numerous as-grown crystals have appeared and were naturally shaped as flat paral- 
lelograms of transparent orange-brown color. The typical sample size was 3 — 5 mm along the longest side of the 
parallelogram. The crystallographic symmetry was checked by X-ray diffraction and found to be triclinic with lattice 
parameters and angles in agreement with the literature^. The plane of the as-grown crystals was perpendicular to the 
a* axis of the reciprocal lattice. The long and short side of the parallelogram was parallel to the c and the (b -I- c) 
direction, respectively. This natural shape of the crystals allowed an easy orientation of the crystals in the magnetic 
field H with H||a*, H 1 {a*,c} or H ± {a*,{h + c)}. 

Additionally, a larger sample was grown from the seed by a temperature gradient method. This sample has a well 
developed plane normal to the a* axis and has its longest dimension parallel to c-axis. For reasons related to the 
sample space configuration this sample was used for measurements with H||a* and H _L {a*,c} only. Samples of 
deuterated PHCC used for the specific heat measurements were grown from the seed in a similar way. 

The samples were characterized by static magnetization measurements with a Quantum Design MPMS system. 
The concentration of impurities was estimated from the low-temperature susceptibility assuming g = 2.0 and 5 = 1/2 
and was less than 4 x 10~^ paramagnetic impurities per copper ion. 

The technique of electron spin resonance is a powerful tool for the study of low energy dynamics of a spin-gap magnet. 
It can access inter-triplet transitions and measures directly the differences between triplet sublevels. Characteristic 
oricntational, frequency-field and temperature dependences of the resonance absorption spectra provide information 
on anisotropic interactions, zero field splittings and the (7-factor anisotropy. 

In case of PHCC additional complications and challenges arise from the low (triclinic) symmetry of the crystals. 
Due to this low symmetry, neither the (7-tensor nor the anisotropy tensor are fixed to any of the crystal axes. This 
means that the axes orientation for all anisotropic contributions should be considered as arbitrary and mutually 
independent. To solve this problem we have taken a series of absorption spectra at different temperatures and at 
different orientations of the applied magnetic field. 
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FIG. 1. (color online) Left panel: Temperature evolution of X-band ESR absorption derivative (arb. units, / = 9.4 GHz, 
orientation close to H ± {a*, c}). Dashed lines are y-magnified by a factor of 10 with respect to the raw data. Low-field part 
of low-temperature absorption curves is y-magnified by a factor of 100 with respect to the raw data. Right panel: Temperature 
dependence of the static susceptibility (filled circles, denoted as "SQUID") and integrated intensity of ESR absorption (filled 
squares, denoted as "X-band", from the analysis of absorption at the frequency of 9.4 GHz; other symbols, denoted as "38 GHz", 
from the analysis of absorption in different orientations at frequencies close to 38 GHz). 

The ESR measurements were carried out in the frequency range from 9GHz to 120GHz and at temperatures 
down to 0.4 K. X-band (9.4 GHz) measurements above T = 2.5 K were done with a Bruker EleXsys spectrometer 
equipped with a He-flow type cryostat and an automated sample goniometer. For higher frequency measurements a 
set of homemade transmission type ESR spectrometers equipped with 6 — 8 T cryomagnets and microwave oscillators 
covering the frequency range up to 120 GHz were used. Measurements below 1 K were performed with a home made 
ESR spectrometer equipped with '^He-pumping cryostat and a 12 T cryomagnet. 

Low-temperature specific heat measurements in an applied magnetic field up to 14 T were performed with a 
Quantum Design PPMS system equipped with a dilution refrigerator. 



III. EXPERIMENTAL RESULTS 
A. Temperature evolution of ESR spectra 

The temperature evolution of the ESR absorption is illustrated in Figures [1] and [51 At high temperatures we observe 
a single-component ESR absorption with a g-factor close to 2, which is typical for Cu^"*" ions. Below approximately 
10 — 15 K the intensity of the resonance quickly decreases. As expected for a spin-gap magnet this decrease can be 
scaled with the decrease of the static susceptibility (see right panel of Figure [T|) . 

As the intensity of the absorption decreases with cooling the absorption line broadens and, around 5 K, it splits 
into several components. This splitting is anisotropic (see right panel of Figure One can distinguish several 
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FIG. 2. (color online) Left panel: Temperature dependence of ESR absorption (/ — 38.04 GHz, H|[a*). Right panel: ESR 
absorption spectra for different orientations at 3.0 K. 

components that continue to fade out with cooling and one component, located close to a g = 2.0 resonance field, 
with an increasing intensity upon cooling. The fading components can be ascribed to a resonance of triplet excitations, 
while an irregularly shaped component with increasing intensity corresponds to defects and impurities. 

The magnitude of the splitting is temperature dependent but appears almost constant for temperatures lower or 
equal to 3 K (sec Figures [1] and . Because of the rapid decrease of intensity of the split spectral features below 3 K 
we performed all measurements for the determination of the splitting parameters at 3 K. 

B. Frequency-field diagrams 

Above the splitting temperature where the ESR absorption consists of a single component the measurements at 
different frequencies proved that the resonance field can be described by an effective g-factor: for H||a* g = 2.05±0.01, 
for H ± {a*,c} g = 2.206 ± 0.002 and for H _L {a*,(b + c)} g = 2.236 ±0.002. 

Frequency-field diagrams for the resonance fields measured at 3. OK are shown on the Figure [31 For all three chosen 
orientations the two main absorption components are located to the left and to the right of the high-temperature 
resonance position. Another, weaker component located at approximately half of the first resonance field (see also 
right panel of Figured]) was observed in the measurements with the larger sample for H||a* and H _L {a*,c}. The 
slope of the main components is determined by corresponding g-factor value. The splitting of these components is 
orientation dependent and is equal to 9.2 GHz, 5.5 GHz and 12.8 GHz for H||a*, H ± {a*,c} and H ± {a*, (b + c)}, 
correspondingly. 

C. Orientational dependences 

The angular dependences of the ESR absorption were accurately measured at X-band frequency (9.4 GHz) for 
T = 3.0 K (below the splitting temperature) and for T = 25 K (above the splitting temperature) where the lincwidth 
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FIG. 3. (color online) Frequency field diagrams for different orientations of the applied field at T = 3.0 K. Open circles: 
high frequency measurements at fixed orientations. Closed squares: resonance field positions for corresponding orientations 
extracted from X-band (9.4GHz) orientational dependences. Curves: perturbative model calculations. 

shows a minimum. Three non-orthogonal rotation axes were chosen for different experiments, namely rotations 
around the a* reciprocal space axis, around the c axis and around the (b + c) direction. The samples were fixed 
with paraffine inside suprasil glass tubes, with the chosen rotation axis parallel to the tube axis. The accuracy of this 
sample mounting technique was limited with respect to the rotations around the tube axis and, thus, the unavoidable 
angle shift was later taken as a fit parameter. 

The positions of the resonance fields shown in Figure S] were determined from fitting the spectra by a single or a 
multi-component absorption, assuming a Lorentzian shape for each of the components. For the data at 3.0 K we did 
not include the weak low-field component (see the line below 0.2 T in the left panel of Fig. [T|), since this line is broad 
and distorted in shape, which makes a precise determination of the resonance field unreliable. 



D. Low temperature ESR measurements 

To check for possible singlet-triplet transitions and antiferromagnetic resonance in the field-induced ordered phase 
we performed additional measurements at the temperatures down to 450 mK using the large sample with H||a*. At 
frequencies of 27 GHz and 33 GHz and for fields up to 12T we did not observe any other resonance absorption at 
450 mK besides that due to defects. On heating above 1.5 K a split ESR absorption, similar to that described in the 
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sample rotation angle (degrees) 

FIG. 4. (color online) Orientational dependences of resonance field in X-band (9.4 GHz) experiments. Selected field directions 
with respect to the crystal axes are shown as determined from the fit. Circles, squares: experimental data at 25 K and 3.0 K, 
respectively. Lines: best fit by the perturbative model described in the text. 

above subsections, was observed. 

E. Critical field measurements 

In order to determine the critical fields at different orientations the field dependence of the specific heat was measured 
at 500 mK (Figure [5]). These measurements were done on deuterated samples of PHCC while non-deuterated samples 
were used for ESR experiments. A sharp maximum of the specific heat marks the transition to the field-induced 
ordered phase. The determined values of the critical fields are: for H||a* fioHc = 8.29 ± 0.05T, for H _L ( — 1, f,0) 
^iqHc = 7.90 ± 0.07T. 

IV. DISCUSSION 

This section is built as follows: first, we will qualitatively describe our data and will show that our observations 
reflect the effects of anisotropic interactions; second, we will make a quantitative analysis of the ESR data and 
will estimate the strength of these anisotropic interactions; flnally, we will discuss the relations between the found 
anisotropic properties at low field with the known anisotropic properties at the critical field. 

A. Qualitative description 

The observed evolution of the ESR absorption spectra is typical for a spin-gap magnet with triplet levels split by 
an effective crystal field. A similar behaviour was observed earlier for other spin-gap mag nets, e.g. TlCuCl;?! and 
NTENB^. 

The decrease of the absorption intensity with temperature corresponds to the decrease of the population numbers 
of gapped excitations. At high temperature the population numbers arc high and excited quasiparticlcs interact with 
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FIG. 5. (color online) Field dependences of the specific heat of deuterated PHCC at different orientations and at T = 0.5 K. 
Curves are guides to the eye. 



each other switching on the exchange narrowing mechanism. This causes a single resonance line observed at high 
temperature. At low temperatures excitations are rare and can be considered as a gas of noninteracting quasiparticles. 
Hence, their spectrum can be affected by anisotropic interactions, which can be described as an effective crystal field 
acting on 5* = 1 quasiparticles. These interactions lift the degeneracy of the triplet sublevels at H = 0, causing 
splitting of the thermoactivated ESR absorption spectrum into several components. This effect is similar to the well 
known splitting of energy levels of an 5 = f ion in a crystal^. The line broadening marks a crossover from the 
high-temperature, exchange-narrowed limit to the low-temperature single-particle limit. 

Besides the thermoactivated ESR absorption which is caused by the resonance transitions between triplet sublevels 
an observation of singlet-triplet transitions or an antiferromagnetic resonance above the critical field could be possible. 
Both of these transitions were observed in the above-mentioned spin-gap systems TlCuCla and NTENP. However, the 
intensity of the dipolar singlet-triplet transition should be zero in the exchange approximation. Thus the observation 
of this transition should require in particular a symmetry-breaking of the ground state by anisotropic interactions. 
The absence of the singlet-triplet transition in our experiment shows that the coupling of these states is small. The 
absence of the resonance absorption in the field induced ordered phase even up to I2T at 450mK (at H = lOT 
T/v w 2.5K—) provides 30GHz (or 0.15meV) as an upper estimate of the gap. 



B. Application of a perturbative model for the description of ESR data 

Different models are used for the description of the low energy dynamics of a spin gap magnet in magnetic field. A 
fermionic field-theory model was proposed for ID systems by TsvelilsiS. Alternatively, a bosonic field-theory model 
was proposed by Affleck^ and later independently developed as a macroscopic model by Farutin^. Particular systems 
can be treated by microscopic approaches^. However, only close to the critical field the differences of all these models 
become important. Since we can not observe singlet-triplet transitions on PHCC, we can not look for these fine 
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FIG. 6. (color online) Orientation of the cartesian basis used for calculations with respect to the crystal axes and orientations 
of anisotropy and g-tensor axes found from the modeling. 

differences and can not reliably judge on the applicability of different models from the ESR data alone. On the other 
hand, since our data correspond to the low-field limit (the critical field is around 8 T, while all X-band resonance 
modes are below 1 T and most of the higher frquency data are below 4 T) we can apply a perturbative approach for 
the description of the observed ESR data. It will yield the values of zero-field splitting of the energy levels and the 
orientations of all relevant anisotropy axes which can be used, if necessary, to obtain parameters of other models. 

This approach considers triplet excitations for a non-interacting S ~ 1 particle moving in a stationary effective 
crystal field. The effective spin Hamiltonian for this particle is: 

n^A + B.-g-S + DSzJ + E{SxJ-SYj) (1) 
where A is an energy gap, 5 is a 5-tensor, D and E are the effective anisotropy parameters and X^i, Y^, Za are 
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FIG. 7. Left: Fragment of the PHCC crystallographic structure (from Ref0)- Vertical direction corresponds to the crystallo- 
graphic c axis. Right: Orientations of the distorted chlorine pyramids. Positions of copper and apical chlorine ions are shown. 
Visible bases of the pyramids are shadowed in grey. 



the anisotropy axes. Due to the triclinic crystal symmetry the orientation of the anisotropy axes X^i, ^Ya, '^a is 
arbitrary with respect to the crystal. Moreover, the 5-tensor axes Xg, Yg, Zg are also arbitrary with respect both 
to the crystal and to the anisotropy axes. We neglect here the k-dependence of the Hamiltonian parameters D and 
E, since at low temperatures only the bottom of the spectrum is populated, while at high temperatures the exchange 
narrowing cancels out the effects of these terms. 

The perturbative model can be solved for eigenenergies for any given magnetic field. The differences of the eigenen- 
ergies will correspond to the resonance frequencies which can then be compared with the experiment. We used a 
standard least-squares routine to minimize the model deviation from the experimental data. 

Angular dependences of the X-band resonance fields at 3 K and 25 K were used for the determination of the 
Hamiltonian parameters. Additionally, the values of the (/-factor and the low-temperature splittings measured in the 
high-frequency experiments with well-defined sample orientation were used as anchors to correct for the angle shifts 
due to the sample mounting. The basis X, Y, Z with Z||a* and X||(— c) was selected to describe the orientation of 
the g-tensor and the anisotropy axes with respect to the crystal (Figure [5]). 

This problem has 14 parameters: the anisotropy constants D and E, the Euler angles of the anisotropy axes a a, 
Pa and 7^, three main values of the g-tensor, the Euler angles of the g-tensor axes, and, finally, three uncontrolled 
shift angles £1,2,3 in the X-band rotation experiments. It turned out, that this general fit yields two of the 5-tensor 
components very close to each other (best fit g-factor values 2.272, 2.062 and 2.039). Consequently, the fit became 
unstable and could not well distinguish the orientations of the g-tensor differing by rotation around high-5- value axis. 
Thus, we reduced the number of model parameters to 12 by assuming an axial g-tensor that can be described by 
two main values and by two polar angles of the (7-tensor axis. The stability of the best fit solution was confirmed by 
running numerous fitting procedures with randomly varied initial approximations of the model parameters. 

The best fit parameters are shown in the Table U and the corresponding frequency-field and angular dependences 
are shown in the Figures |31 |H The orientation of the anisotropy axes and g-tensor axis are illustrated at the Figure 
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TABLE I. Best fit parameters in the perturbative model. 



anisotropy 


D, MHz 


— 7900 ± 280 


constants 


E, MHz 


-1340 ± 190 


g-tensor 


511 


2.280 ± 0.015 


main values 




2.048 ± 0.007 


Euler angles 


OLA 


(-43.5 ± 1.8)° 


of anisotropy 


Pa 


(—105.6 ± 1.5) 


axes 




(—12.5 ± 7)° 


Polar angles of 




(105 ±3)° 


g-tensor axes 




(54.9 ± 1.5)° 


shift 


ei 


(-10.9 ±3.5)° 


angles 


£2 


(41 ±2.5)° 




£3 


(-4.0 ±3.5)° 



El The sign of anisotropy constants can be unambiguously determined by comparing low-temperature intensities of 
the spht absorption components in different orientations (Figure [2]) : a more intense absorption corresponds to the 
transition from the lower sublevel. 

The found axiality of the g-tensor can be explained by the details of the microscopic structure of PHCCj^ Each 
copper ion is surrounded by five chlorine ions, forming a distorted square-based pyramid (Figure [7]). The main axis of 
this pyramid is the same for all copper ions (differently oriented pyramids are connected by inversion). The direction 
of the vector linking the copper to the apical chlorine (position Cl(2) in the denominations of Ref|j) is very close 
to the found direction of the g-tensor axis: its polar angles in the chosen reference frame are Qcu-Ci = 111° and 
ipcu-ci = 57°. 



C. Anisotropic corrections and critical field. 

Our specific heat measurement and earlier measurements of other authors demonstrate the anisotropy of the critical 
field of about 1 T (see Table HI)) . Additionally, neutron scattering experimentsiii^ have determined the orientation of 
the order parameter above the critical field for H||b: the order parameter was found to lie within (7 ±4)° of the (ac) 
plane with an in-plane component within (6 ±4)° of the a-axis. Both of these effects are determined by the anisotropic 
interactions in the spin system and thus provide an independent check for the anisotropy parameters found above. 

The critical field depends strongly on temperature: its values at 65mK and at 2.8K differ by w 2T— . Thus, it is 
necessary to ascertain whether the anisotropy parameters determined at 3. OK can be used in the analysis of the low 
temperature data or not. 

The temperature dependence of the critical field originates from the renormalization of the energy spectrum due 
to the repulsive interaction of triplet excitations. This renormalization depends on the triplet population numbers 
and therefore affects strongly the energy spectrum in the vicinity of the critical field, where the energy gap is small. 
Our ESR measurements, on the contrary, were performed in small fields where the energy gap value is still high and 
the triplet population numbers are small even at 3. OK. This is confirmed by the experimental observation that the 
observed splitting of the resonance absorption does not change significantly at further cooling below 3K. Therefore, 
the effect of triplet repulsion does not affect the values of anisotropy parameters found in the previous subsection 
and these values can be used for the analysis of the low temperature data for the critical field anisotropy and order 
parameter orientation. 

In the model calculations of this subsection we will use g-factor values and (when necessary) zero-field splittings 
as determined in our ESR experiments. Since the gap value A can not be determined with the similar precision, we 
will tune it to obtain correct value of the critical field in the chosen orientation. We choose the value of the field at 
which the antiferromagnetic Bragg refiection appears^ at H||b (/ioi^c = 7.47±0.01T) as the most reliable anchor for 
this purpose. 

First, there is a simplest model, suggested in Ref. [l|, that treats the anisotropy of the critical field as originating 
from the gi-factor anisotropy only. This model was backed by the observation that the ratios of the critical field to 
the saturation field were found to be the same in different orientations. The values of the g-factor in the particular 
orientations of the magnetic field calculated from our model are given in Table [ill One can see that the highest 
critical field is indeed observed in the direction with the smallest g- value (H||a*). The critical field in the arbitrary 
orientation can be expressed via the ^-factor value as He = (5(H||b)i7c(H||b))/g. The calculated values are in 
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TABLE II. Comparison of the measured parameters of the excitation spectrum and calculated parameters in different models. 
Values of g-factor are given for reference as calculated for particular orientations using a best fit perturbative model. For the 
precision of the model estimations see footnote 0. 





Experiment 


g-factor 
anisotropy 
only 


pertur- 
bative 
model 


macro- 
scopic 
model 


Zero-field 
energies, 
GHz 


247^ 


235.6 


233.6 
236.2 
242.8 


234.0 
236.6 
243.2 


/io-ffc,T: 










H± (-1,1,0) 
{g = 2.183) 


7.90 ± 0.07^ 


7.71 


7.74 


7.70 


H||a* 
{g = 2.064) 


8.29 ± 0.05^ 
8.03 ± 0.07^* 


8.15 


8.27 


8.19 


H||b 
{g = 2.253) 


7.47 ± 0.01^ 
7.58 ±0.03^ 
7.40 ± 0.08 


7.47 


7.47 


7.47 


Hl|c 
(S = 2.122) 


7.6 ± 0.2_^ 


7.93 


7.96 


7.89 



Estimated errors of the model calculations: Absolute value of zero field energy it2GHz. Difli'erence of sublevels energies ±0.4GHz. 
Critical fields; itO.OST. This error is mostly due to the g-factor uncertainty. 
Ref(3, deuterated PHCC, inelastic neutron scattering, T=60mK 
present work, deuterated PHCC, C(H) at T=0.5K 
Ref U magnetization measurements at 0.46K 
" Ref 3, deuterated PHCC, elastic neutron scattering, T=65mK 



reasonable agreement with the experiment (see Table HII . 

However, this simple model is not exact, since we observe zero-field splitting of energy levels. The effect of zero- 
field splitting on the anisotropy of the critical field can be roughly estimated from the value of the main anisotropy 
constant as ^qAHa ^ hD/{gfj,B) ~ 0.3T. This has to be compared with the effects of the g-tensor anisotropy: 
fioAHg ^ {Ag/ g)fioHc ~ 0.9T. These values are similar and, especially since the axes of anisotropy and the axes of 
g-tensor do not coincide, it is not possible to neglect a priori one of the effects for an arbitrary field direction. 

As a second test model, we extrapolate the results of the perturbative model up to the critical field. This ex- 
trapolation can be justified only for one-dimensional systems^. However, this model can be considered as a simple 
approximation taking into account both the effects of g-tensor anisotropy and zero-field splitting of triplet sublevels. 
The results of this extrapolation are also collected in Table [Hi 

Finally, the problem of the triplet level field dependence can be treated more strictly within the frameworks of 
a macroscopic (bosonic) approach^. This model treats the excitations of a spin-gap magnet as oscillations of the 
vector field r) with the Lagrangian 



C=l{ri + l[v^VL]f~lAv' + Crei (2) 

where A is the exchange constant describing the energy gap, 7 is the gyromagnetic ratio and Crei is a term containing 
anisotropic relativistic corrections. The leading corrections terms at low fields are the effective anisotropy terms, 
responsible for the zero field splitting, which can be written in the reference frame of anisotropy axes {Xa,Ya, Z a): 



Crell = \hl{vl_, + - 27]|J + ^^2(??L " ??f-J (3) 

where 61^2 are the effective anisotropy constants. These terms allow to describe the zero field splitting, see e.g. Refl^. 

To describe additional field effects we have to include here next order terms. The low symmetry of the crystal allows 
numerous relativistic and exchange-relativistic invariants such as fjafjp, [ry x rf\aHp, {iiri)Harif}, H^rjarjfj, rfHaHp 
etc. I.e., in general we have to consider a veritable zoo of possible terms which overburdens the problem beyond any 
hope to produce a compact solution. 
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To simplify this problem we will assume that the main field effects arise from a single-ion g-tensor anisotropy. 
Thus, due to an almost tetragonal local symmetry, these effects can be considered as having an axial symmetry in the 
(7-tensor reference frame. Additionally, we will require that for the model test case of 61 = 62 = these terms should 
provide linear E(H) dependences for all three modes, which is natural for the microscopic model with anisotropic 
single-ion g-tensor only. 

This leads to the following invariants combination: 

Lr,n = ai^Hl^rf - j^Hri)Hz^vz, + iHz, [il x 77] z,) (4) 

where the parameter ^ = {g\\/gi_ — 1) reflects the g-factor anisotropy. Up to higher order corrections, equation ^ 
is equal to the replacement of the problem with an anisotropic g-tensor g and a certain real field direction H by the 
formally equivalent problem with the isotropic g-factor 50 and a transformed "effective" field He// = {l/gQ){gii). 
Three terms in Eqn.(j3]) are responsible for different effects and none can be omitted. 

The dynamics equations can be obtained from the final Lagrangian C + Creii + Crei2 using a standard variational 
technique with special attention to the frames of reference (see Appendix). The solution of these equations yields the 
field dependences of the energy levels, whose differences would yield transition frequencies. The parameters of the 
macroscopic model can be directly expressed via perturbative model parameters with only one adjustable parameter 
A: 

A = 4:7r^{A^ + {A + D + Ef + {A + D - Ey)/3 
bi = 4:7r^{2A^ - {A + D + Ef - {A + D - E f)/6 
&2 = 4:7r^{{A + D + Ef - {A + D - Ey)/2 

7 = 27r(gx/2.0)7o 

e - {g\i/gi_ - 1) 

here 70 = 2.80GHz/kOe is a free electron gyromagnetic ratio. The factors of 47r^ or 2tt appears here since the dynamics 
equations (jA.2|) are written for angular frequency uj, while the parameters of the perturbative model were expressed 
in the ordinary frequency units. 

The parameters values corresponding to ^o-ffc(H||b) = 7.47T arc: A = 2.236 • 10'^^ l/scc^, 61 = 4.972 • 10^^ l/scc^, 
b2 = -2.490 • 10^2 l/sec2, 7/(27r) = 2.867GHz/kOc, ^ = 0.1133. 

We have found that the model curves fits our data perfectly and are undistinguishablc from the perturbative curves 
at low fields. The predicted values of critical fields are also gathered in the Table HIl 

The macroscopic model allows a straightforward analysis of the order parameter orientation above the critical field. 
We have to add to the final Lagrangian fourth order exchange term —Brf which will determine the magnitude of the 
order parameter al H > However, as long as all anisotropic terms arc of the same order in rj, the orientation 
of the order parameter docs not depend on its amplitude and, consequently, docs not depend on the value of the 
exchange constant B. Thus, for the sake of determining the order parameter orientation only, the B constant can be 
simply thought of as a mathematical meaning to exclude the divergence of the order parameter above He- Using the 
model parameters described above we have found that the polar angles of order parameter at H||b are Qq.p. = 34° 
and (po.p. = —20°. This corresponds to an angle of 16° with the (ac)-plane and to an angle of 28° to the a-axis. 
The in-plane component of the calculated order parameter forms an angle of 23° with the a-axis. Although there is 
a qualitative agreement, the calculated values differ from those found in the experimenljiii^. 

All of the models provide a reasonable qualitative agreement with the experiment. The macroscopic model seems 
to be physically the most justified, since it describes correctly the zero-field structure of energy levels and can be 
applied up to the critical field and above. Quantitative disagreements with the experiment can be caused by some 
neglected factors: the slight change of anisotropy constants at cooling below 3K or the effects of additional anisotropic 
contributions above He- 

V. CONCLUSIONS. 

A detailed investigation of the low-energy spin dynamics in the spin-gap magnet PHCC revealed the presence of 
anisotropic interactions causing a splitting of the triplet sublevels and leading to a 5-tensor anisotropy. The effects 
of these interactions were described quantitatively within a perturbative approach in small fields. The parameters 
of anisotropy comprising the main g-tensor values, the zero-field splitting of sublevels and the orientations of corre- 
sponding axes with respect to the crystal were found. Additionally, an application of the macroscopic (bosonic) model 
to the case of a spin gap magnet with a relevant 5-factor anisotropy was developed. 
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Appendix: Dynamics equation for macroscopic (bosonic) model with axial g-tensor below He- 

The complete Lagrangian of the problem can be obtained from Eqns.®,®,!!!]): 

= i(77 + 7[^xH])2-lAr,2 + 

+ ^jHz^iv X 'n]z^ (A.l) 

The same coefficient ^ in the last three terms is a consequence of the microscopic model that assumes a single-ion 
g-tensor anisotropy with an axial g-tensor. Due to the low symmetry of the crystal, terms responsible for the effective 
anisotropy (6i and 62 terms) and for the g-factor anisotropy (^-terms) are written in different cartesian bases. When 
varying the Lagrangian over rj these terms have to be recalculated to the same basis. We used the anisotropy axes 
basis for the derivation of the dynamics equation. The resulting dynamics equation on r] is: 



r) + 2-f[r] X H] - -f'^H^T] + 7^H(r7 • H) + At? 





1 






1 


VXa \ 


bi 




VYa 






-VYa 




\ 


-^VZa I 




\ 


/ 



C7'(HZ,)((r,Z,)H + (r7H)Z,) 
2^7(HZ,)[Z, xr,]=0. 



(A.2) 



This equation yields a secular equation for the precession circular frequencies uji^2.3 that describes the field depen- 
dence of the energy levels. The latter is a cubic equation for uj"^ and can be solved for any field direction. 
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